Stick-slip dynamics of an oscillated sessile drop 
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The dynamics of an oscillated sessile drop of incompressible liquid with the focus on the con- 
tact line hysteresis is under theoretical consideration. The solid substrate is subject to transverse 
oscillations, which are assumed small amplitude and high frequency. The dynamic boundary con- 
dition that involves an ambiguous dependence of the contact angle on the contact line velocity is 
applied: the contact line starts to move only when the deviation of the contact angle exceeds a 
certain critical value. As a result, the stick-slip dynamics can be observed. The frequency response 
of surface oscillations on the substrate and at the pole of the drop are analyzed. It is shown that 
novel features such as the emergence of antiresonant frequency bands and nontrivial competition of 
different resonances are caused by contact line hysteresis. 
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PACS numbers: 47.55.D-, 47.35. Pq, 47.55.np, 46.40.-f 
I. INTRODUCTION 

Rapid development of microtechnologies over last 
decades has manifested great interest in theoretical as- 
pects of contact line dynamics. The abihty to predict 
the motion of contact line and hence to control wet- 
ting processes becomes of paramount importance for 
applicationsJ^i^ Despite noticeable progress in under- 
standing the steady motion of the contact line, the un- 
steady motion remains significantly less explored and in- 
volves a number of important open questions. Of special 
interest is the role of contact angle hysteresis, which for 
unsteady motion may become a crucial feature in ob- 
taining the proper picture of the contact line motion. In 
the present study, we address this issue in the context of 
oscillated sessile drop. 

The dynamics of a drop on oscillated substrate has 
been considered for many years, see recent surveys in 
Refs. j^.^. Recent experimental studies have shown that 
the contact line hysteresis can lead to such nontrivial ef- 
fects as stick-slip dynamics of the contact line^, climb- 
ing motion over inclined substrate^ and motion over 
gradient^i^ and thermal gradient^" surfaces. These ex- 
perimental observations raise a natural question about 
the role of the contact angle hysteresis, which is easy to 
pose but rather difficult to answer. In most theoretical 
studies, the contact angle hysteresis has been either com- 
pletely neglected^iiiii^ii^ or treated in an oversimplified 
way, where the drop is similar to an oscillator with solid 
frictioujASiiiii^ii^ii^ Although these solid-friction models 
reflect the qualitative picture of the stick-slip process, 
they do not provide satisfactory understanding of the 



phenomenon. 

To obtain deeper insight into the physics of the 
stick-slip motion, a boundary condition suggested by 
L. M. Hocking— can be applied. This condition, which 
captures principal features of the contact line motion, in- 
volves an ambiguous dependence of the contact angle on 
the contact line velocity: 



9C 
dt 



A(x - Xo), X > Xo, 
0, Ixl < Xo, 

A(x + Xo), X < -Xo- 



(1) 
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Here, functions C and x describe the deviations of the 
free surface and the contact angle from those in equilib- 
rium, respectively, and xo is the critical value defining 
the contact angle hysteresis (see, e.g.. Fig. [T]). The fac- 
tor A, which has the dimension of velocity, characterizes 
interaction between the substrate and the liquid and is 
referred to as the wetting or the Hocking coefficient. 

The particular case of xo = 0, in which the con- 
tact line velocity d(/dt oc x, describes no contact an- 
gle hysteresis. Different practically important situa- 
tions can be addressed by changing A. In terms of the a 
corresponding dimensionless parameter, for instance. A, 
as in relation ([S]), these situations range from the com- 
pletely pinned contact line, A (the contact angle can 
change) to the opposite case of the fixed contact angle, 
A — > oo (the contact line is freely moving). 

In the present study, we consider the dynamics of 
a hemispherical drop on a normally oscillated solid 
substrate. We apply condition ([T]) without compromise 
and focus on the role of the contact line hysteresis, which 
has been recently measured in a similar setup. ^ Based 
on this approach, we are not only able to quantitatively 
describe the stick-slip process but particularly to reveal 
a new interesting feature in the contact angle evolution. 
This finding is very much reminiscent of the experimen- 
tal observations, which cannot be explained in terms 



2 



2 ; 


1 ambient 


A ^ 

(f drop 

T '\ 




1 substrate | 



FIG. 1: Problem geometry. A hemispherical drop on the 
transversally oscillated substrate. 



of the previously suggested theoretical models. The 
paper is outlined as follows. We start with the problem 
statement in Sec.|lll Section Hill provides the description 
of the method we use to treat the problem. The obtained 
results are discussed in Sec.lIVIand summarized in Sec. El 



II. PROBLEM STATEMENT 

Consider a sessile drop of incompressible liquid of den- 
sity p and kinematic viscosity v, see Fig. [1] We are in- 
terested in the situation of a gaseous ambient, where its 
density is much smaller than that of the drop and there- 
fore can be neglected. We assume that the solid substrate 
is subject to transverse oscillations with an amplitude a 
and a frequency uj. We admit that the drop is enough 
small so that its shape is hardly distorted by gravity. 
This assumption ensures that the equilibrium drop sur- 
face is hemispherical to very high accuracy, with radius 
R, and the equilibrium contact angle equals 7r/2. 

The amplitude of external driving is considered small 
in the sense e = a/ R <C 1 and the frequency of the sub- 
strate oscillations is high enough: ujR^ ^ v. At such 
frequencies, viscous boundary layers, which arise near 
the rigid plate and near the free surface, become very 
thin. In other words, the frequency restriction allows us 
to neglect viscous dissipation in the liquid, which ensures 
that the approximation of inviscid liquid is justified. ^'^ On 
the other hand, the frequency lo is assumed comparable 
with the eigenfrequencies a;„ of shape oscillations for a 
spherical drop of radius R. For our consideration, only 
even eigenfrequencies are of interest, which spectrum is 
defined by the relation = 2n(2n - l)(2n 2)cr/(pi?3), 
where a is the surface tension. 

Because of symmetry, we use the spherical reference 
frame with the coordinates r, z?, a and the origin at the 
center of drop and restrict our analysis by the axisymmet- 
ric problem. In the accepted approximations, the fluid 
motion is irrotational, which makes it convenient to intro- 
duce the velocity potential ip. As a result, the dynamics 
of the liquid is described by the Bernoulli equation and 
the incompressibility condition. Let r = R + t) be 
the instantaneous locus of the distorted free surface, see 
Fig. [H To make a comparison with the nonhysteresis 



study^ simpler, we measure the time t, length r, velocity 
potential ip, the deviation of the pressure field p from its 
equilibrium value and the surface deviation ( in the scales 
of -\/ pR^/a, R, a^/aJpR, aa/ R^ , and a, respectively. As 
a result, the dimensionless boundary value problem is 
defined by (intermediate steps can be found in Ref. 0) 



at 

■K dC { '^('^ "'^o), 7 > 70, 

r = l,d^-: { 0, 

I A(7-|-7o), 7 < -70- 



(2) 
(3) 
(4) 

I7I < 70, (5) 



Here, the differential operator 
sin "d dd 



sini9 



5C 



and 7 — ~{dC,/ dd)\^^T^/2 is the dimensionless deviation 
of the contact angle from its equilibrium value, which for 
the sake of brevity will be called the contact angle. 

Boundary condition Q ensures impermeability of the 
substrate for the liquid. The kinematic and the dynamics 
conditions at the free surface are presented by Eq. 
In contrast to the previous work^ we impose a more gen- 
eral Hocking condition at the line of contact of the three 
phases, as given by Eq. Thus, boundary value prob- 
lem ([2])- ([5]) is similar to the previous study, except for 
this hysteretic condition. 

The problem involves three dimensionless parameters 



puj'^R^ 



A = X\ 



70 = 



xoR 



(6) 



which have the meaning of the squared external fre- 
quency rescaled with respect to the eigenfrequencies ujn, 
the wetting (or Hocking) parameter, and the critical 
value of contact angle, respectively. 

As the frequencies lo and ujn have been assumed com- 
parable, the parameter f2 is finite. A similar argument 
refers to the parameter 79. We focus on the case of well 
polished substrate, which implies that the threshold value 
Xo is small. Being the ratio of two small parameters, xo 
and e, the parameter 70 is treated as finite. This observa- 
tion indicates that the contact-line hysteresis is expected 
to be non-negligible even for small-amplitude oscillations. 

We emphasize that condition ^ admits a number of 
particular cases. In the limiting case of perfectly polished 
surface, 70 — > 0, the boundary condition ^ is reduced 
to its simplified modificationi^ Particularly, the contact 
angle remains fixed as A — > 00, whereas the opposite case, 
A — > 0, describes the contact line pinned. We note that as 
in the latter case, the same dynamics refers to the limit 
of large values of 70. 
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Finally, we stress the nontriviality of our considera- 
tion. Although the amplitude of oscillations is considered 
small, e <C 1, which allows us to linearize the govern- 
ing equations and simplify the boundary conditions, the 
overall problem is nonlinear. The nonlinearity is ensured 
by posing the Hocking condition ((5]) and comes into play 
through the parameter 70 ~ a^^, therefore the solution 
is eventually amplitude dependent. 



III. METHOD OF SOLUTION 

To treat the formulated problem, we note that the ve- 
locity potential satisfies Laplace's equation ([2]) and there- 
fore can be presented as a series in the Legendre polyno- 
mials. In view of impermeability condition ([3]), only even 
terms are nonvaninshing in this expansion. We retain the 
terms regular at the origin and present the solutions for 
the velocity potential and the consistent solutions for the 
surface deviation and the pressure in the form 



P 



n=l 



Cn{t)P2n{0y 

2n 



(7) 
(8) 



Po[t) — y — it zcosiU, (9) 

Tl=l 



2n 



where 9 = cosd. Note that the summation for ( starts 
from n ~ I. The term with n = is set to zero to ensure 
that the drop volume is conserved. The zeroth harmonics 
ifoit) and po{t) = —{d(po/dt) are nonvanishing spatially 
independent functions, which describe spatially uniform 
pulsations of the velocity potential and the pressure, re- 
spectively. The function po is determined by the require- 
ment of conservation of the drop volume, the term ipo is 
unimportant for the further analysis. 

To distinguish between the time intervals of the contact 
line in motion from those in standstill, we next stick to 
the following notation. Let t = to he the time moment 
when the contact line stops to move and t = ti he the 
switching time when it proceeds with the motion again. 
We next obtain the solutions describing the dynamics 
at these time intervals separately and then show how to 
match the solutions. 



A. Pinned contact line 

We now consider the time intervals characterized by 
small values of the contact angle, [7! < 70. As it follows 
from the Hocking condition ([5]) , the contact line remains 
fixed during this phase of evolution (A = 0) and we can 
use the solution obtained before.— The only point we have 
to care about is that in our situation the contact line is 
fixed not necessarily at r = 1 but at a slightly differ- 
ent position characterized by C = C/ 7^ Oi which is easy 



to take into account. As a result, the solution for ( is 
presented as a superposition of the eigenmodes Co'"\ ^ 
particular solution caused by the external force, ^p, and 
the term cx Q, which corrects the contact line position, 
respectively: 



771—1 

+C/(l-20). 



(t-to) 



(10) 



The form of the last term is chosen such that the overall 
expression for ^ ensures the constant volume of the drop. 
A similar ansatz for the velocity potential reads 



m— 1 

+ipp{r,e,t). 



(11) 



Here, Dm are the complex amplitudes of the eigenoscilla- 
tions to be determined as described in IHI CI and ojm are 
the eigenfrequencies for the drop with the pinned contact 
line. These eigenfrequencies are defined as the roots of 
the transcendental equation 



/(c^,0)=0 



with the function 



n=l 



where 



^ (4n+l)P2n(0) 

(2n~ l)(2n-f 2)' 
nl = 2n(2n- l)(2n + 2). 



(12) 



(13) 



(14) 
(15) 



The values r2„ have the meaning of the dimensionless 
eigenfrequencies of even eigenmodes for the spherical 
drop oscillations. 

The eigenfunctions for the problem with the fixed con- 
tact line (A = 0) are known to be^ 

00 

C^"'(0) = 2 E nAmnP2n{e) = Bmf{uJm,0), (16a) 

00 

4™^(r,0) = iu:m AmnP2n{ey (16b) 

n=0 

with the coefficients 

= „>0 ,17) 



E 



«»»^P2n(0) 



(18) 
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Here, we introduce the normalization condition and point 
out the orthogonahty of the eigenfunctions 



4k) 



(19) 



with (5=1 for m = k and 6 = otherwise. 

For the problem of forced oscillations the solutions can 
be expressed as 



Cp{9,t) = Rc[Cp(%*"* 
ipp{r,9,t) = Re [ipp{r,9)e 
with the complex amplitudes 



:atl 



E 

.n=0 



(2n- l)(2n + 2) 



(20a) 
(20b) 



(21a) 
(21b) 



where 



g{e) = Fe--[i-e\n{i 



En 

F 



1 + (2n - l)(2n + 2)F 



imo)5i 



2na„F2«(0) 



(22) 
n > 0,(23) 

(24) 



We indicate that sum included in relations 
(|16ap and converges very slowly. From the compu- 
tational point of view, its evaluation can be significantly 
improved if a more suitable form is used. By taking into 
account the expansion 



(25) 



n=0 



sum (|13p is presented in an alternative way 
1 



2 O^nP2ni0) 

2 + ^ - ^ 

n=l " 



(26) 



which compared with the original representation ([Te 
provides much faster convergence. Note that a similar 
procedure can be applied to the sum in relation (jlSp . 



B. Moving contact line 

We next deal with the time intervals of supercritical 
values of the contact angle, I7I > 70, when the contact 
line is no longer fixed. We might build the solution in the 
form of series as we did in Sec. IIII Al This way is how- 
ever not worth implementing because the corresponding 



eigenvalue problem is not hermitian and hence no or- 
thogonality condition as in exists. As a result, this 
approach becomes computationally inefficient and pro- 
vides no advantages any more. What we do instead, is 
addressing Eq. ([2]) for the pressure, taken at r = 1, which 
allows us to figure out how the contact line evolves. By 
using the dynamic boundary condition, Eq. we ex- 
clude the pressure and apply ansatz ([7]). As a result, we 
obtain the inhomogeneous Legendre equation 



d_ 



1- 



dc 
de 



Cn{t)P2n{e) 



2n 



(27) 



The solution of this equation is given by 



C{e,t) = -r^^cosm^ 



anP2n{0) 



-^^ (271- l)(2n-t-2) 

f> Cn{t)P2n{e) ( i\ 

- E2.(2.-l)(2n + 2)+n'"2j '^^^^ 

where the term cx 7 is the general solution of the ho- 
mogeneous equation and the first two terms present a 
partial solution of the inhomogeneous equation. We note 
that the integration "constant," 7 = 7(i), satisfies the 
definition of the contact angle 



961 



7- 



(29) 



For this reason, the time-dependent function 7(f) is re- 
ferred to as the contact angle. 

By making comparison of expressions ([7]) and (|28p and 
using relations (|15p and (|25|) , we derive a set of ordinary 
differential equations for the expansion coefficients Cn 



Cn + ^1,Cn = ^^^an7 — 2nil^ Un COS fit, 



(30) 



which are coupled to each other through "f{t). As the 
function j(t) is unknown, an additional relation is re- 
quired to close the system. The formulation of the prob- 
lem is completed by rewriting the Hocking condition , 
which yields 



7(t) 



Sit) 
S{t) 



-70, 7 > 70: 

70, 7 < -70, 



with the auxiliary function 

00 

S{t) = - ^C„(t)P2„(0). 



(31) 



(32) 



Thus, at the stage of evolution with supercritical con- 
tact angles, I7I > 70, we numerically solve the system of 
ordinary inhomogeneous differential equations pop to- 
gether with algebraic coupling relation PT|) . 
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C. Matching the solutions 

To obtain the solution of the full problem, we have 
to match the solutions obtained in Sees. IIII Al and IIIIBI 
The regime with the motionless contact line, which is 
characterized by subcritical values of the contact angle, 
I7I < 7o, is described by expressions ^TU\i and (fTT|) with 
the unknown complex- valued coefficients Dm- At super- 
critical values of the contact angle, I7I > 70, the con- 
tact line keeps moving. At this time interval, we treat 
Eqs. (l30 |) -(|32 |) numerically, in terms of functions C„(i) 
and Cn{t)- Next we provide relations between the coef- 
ficients Dm and C„, C„ valid aX t = and < = ti, at 
which the regimes are switched. At these moments, the 
contact angle reaches its critical value, 7 = 70, and the 
two solutions coincide. 

Consider first the moment t — Iq, when the switchover 
from the regime with the moving contact line to the one 
with the pinned contact line occurs. Given the values 
Cn(to) and C„(to) obtained at the previous phase of mo- 
tion, we have to determine coefficients Dm ■ We multiply 
Eqs. Uni) and mi) by </?^™^(l, e) and Co"^(6'), respectively, 
take into consideration integral condition (|19p . and ob- 
tain the real, Dm"^ , and imaginary, Dm , parts of Dm 



A 



^™ - V-^(C„-27ii?„cosmo) 



71=1 



1 ^ 

LUm ^ 4n + 1 V 

Tl — 1 



Here, the values LOm are determined as the roots of 
Eq. p2|) and the coefficients Amn , Bm , and £"„ are given 
by relations (fT7|) . (fT5)) . and P^l) . respectively. 

This transformation identifies solutions pil)) and (fTTj) 
uniquely and allows us to determine the moment of the 
backward switchover, t — ti, when the contact line starts 
to move again. To find out this moment, we numerically 
solve the algebraic equation with respect to ti 



l{ti) - 70 = 0. 



(33) 



To evaluate 7 we use Eq. (l29l) with the solution for C, 
Eq. (flOl) . where the sum in expression (|16a|) should be 
taken in the form (pS)) . As a result, we obtain 



C„(<i) and Cn{ti)- We equate solutions and (fTTj) to 
those in relations ([7]) and ([5]), which are taken at r = 1 
and t = ti. All the terms are presented as series in the 
Legendre polynomials. This can be done with the aid of 
expressions (fT6)) . pO|. and pT|) . The coefficients on the 
left and right hand sides must be equal, which yields 



d^m'{ti)Amn + E,, COS 



m— 1 

00 



UJmd^rii (^l)^rTm + siu Qti 



Thus, we start from these initial values and solve 
Eqs. (I30l)-(l32l) until the condition I7I < 70 is fulfilled. 
Suppose that 7 crosses the critical value 70 between the 
time steps k and (fc -I- 1). We need to find the moment 
t — to with the accuracy consistent with the integration 
scheme and evaluate the values C„(to) and Cn{to)- To 
estimate to we implement an elegant method suggested 
by Henoni^ The idea is to introduce a characteristics 
that changes its sign while 7 is crossing the value 70. A 
suitable quantity satisfying this requirement can be the 
function S{t), see relation (15^ . which turns to zero at 
7(^0) = 7o- Thus, ^0 as well as the values C„(io) and 
Cn{to) are obtained by making one additional corrective 
integration step S{tk) from the values C„(ifc) and Cn{tk), 
which are available at time step k. 

To have this idea implemented, we proceed from 
Eqs. (|5D|) to the system of differential equations with the 
new independent variable S and the time t = t{S) as an 
additional variable 



~dS ~ ll' 

dCn {oinl ~ C'„) — 2nfl'^an COS rii 



dS 

dt _ 1 
dS ~ H' 



H 



Here, we have introduced the Henon function 



dS 1 " , , 



(35a) 
, (35b) 
(35c) 

(36) 



n=l 



-f{t) = n^F cos nt - 2Cf + Y Bmd^^ (t). (34) 

m— 1 

Here, we have introduced a complex valued function of 
time dm{t) = DmQxp[iu!m{t — to)] with the real and 
imaginary parts denoted as dm'^ and dm , respectively. 
Thus, having obtained the value ti, we are ready to 
proceed to the next situation. 

We now turn to the consideration of the moment t — ti, 
when the motionless contact line starts to move. Before 
we treat Eqs. ([50)) - (|5^ . we need to evaluate initial values 



where C„ can be expressed from Eqs. ((30)) . 

Finally, we integrate Eqs. ([35]) along with relation ([3T|) 
until S changes its sign. We note that while making the 
regular integration steps one sets H = 1. The correc- 
tive integration with the step S{tk) is made with H in 
the form ([55]) . which after the correction corresponds to 
5 = or equivalently to t = ti and hence provides the 
required C„(ii) and C„(ti). 

Thus, we started with the consideration of the moment 
t ~ to, provided the way of proceeding to the moment 
t = ti and then to the next moment t = to- To this end, 
we have obtained the solution over half of the period and 
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the described matching procedure can be successively re- 
peated to obtain the solution at longer times. 

In numerical calculations, infinite number of eigen- 
modes in Eqs. pO|) and (fTTjl was truncated to retain M 
terms. The presented results were calculated for M = 10. 
The control computations with the number of the eigen- 
modes with up to M = 20 have indicated no change in 
the results. However, a further increase of M leads to 
the emergence of unphysical oscillations. The number 
A'' of Legendre harmonics retained in Eqs. ([Zl)-® and 
Eqs. (|16p and (|2ip was chosen to be 100 in most of calcu- 
lations. In order to check the accuracy of calculations we 
performed a number of tests with N = 150 and = 200, 
which gave very close results. 

IV. RESULTS AND DISCUSSION 
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FIG. 2: Characteristics of the steady-state oscillations at A = 
1, 70 — 1, and = 3. Panel (a): Evolutions of the contact 
angle, 7(t), deviations of the free surface at the pole, 
and at the contact line, Co(i)- The horizontal dashed lines are 
the lines 7 = ±70. The filled areas display the time intervals 
of the subcritical contact angles, I7I < 70. Panel (b); The 
Fourier power spectrum evaluated for 7(t). 



We start our discussion by recalling the facti^ that de- 
spite the neglected mechanism of viscous dissipation, the 
Hocking condition itself is dissipative. Exceptions are 
the particular case of the pinned contact line (A — > 0) 
and the contact line freely moving (A — > oo) over the 
perfectly polished (70 0) substrate. Because the sys- 
tem under consideration is generally dissipative, any ini- 
tial state approaches the terminal oscillatory state after 
a certain transient. In other words, any phase trajectory 
is landing at a limit cycle. Because in the case 70 ^ 
and A = 0(1) the decay rate is comparable with the 
frequency of oscillation^ the transient time is estimated 
to be a few periods of oscillations. For these reasons, we 
are mostly interested in the properties of the steady-state 
oscillations. 

To get an impression about the dynamics, we look 
at the steady-state oscillations of the contact angle, see 
Fig. Ufa). When one sees the evolution of 7, it might be 
thought of simple linear oscillations. However, despite a 
relatively simple form of the observed signal, the oscilla- 
tions are nonlinear, which is well seen from the Fourier 
spectrum of the signal. Fig. (^Jb). We note that although 
the driving frequency dominates, a few higher harmonics 
are nonvanishing. Another feature one can readily notice 
is the absence of the even harmonics, which reflects the 
fact that the response is presented by an antisymmetric 
function. This antisymmetry of the steady-state oscilla- 
tions can be seen directly from our mathematical model. 
Indeed, in the terminal state the system is invariant with 
respect to the transformation 

t^t+^, C^-C, (37) 

which is easily understood if one takes into account 
two circumstances. First, the governing equations and 
boundary conditions, Eqs. ©-©J are hnear, if consid- 
ered at intervals with the fixed and moving contact line 
separately. Only the periodic switching between these 
regimes makes the problem nonlinear. We also note that 
the problem, Eqs. (H])-®, is inhomogeneous. However, 



the inhomogeneity cx cos fit changes its sign under time 
transformation (|37p . as required. Second, we take the 
same threshold value 70 used for the advancing and re- 
ceding motion of the contact line. Our numerical tests 
with 70 and 71 ^ 70 for the thresholds of the advancing 
and receding motion, respectively, have confirmed this 
statement. For the distinct threshold values, we detect 
nonvanishing even harmonics, which contributions to the 
power spectrum become more pronounced as 70 and 71 
become more distinct. In the opposite case of 71 — > 70, 
the even harmonics die out and we come back to the 
perfect antisymmetry, as in relation (|37|) . 

Along with the contact angle, we measure the devi- 
ations of the free surface from its equilibrium position 
at the pole of the hemisphere, Cp(i) = — liO and 
on the substrate Co(<) = ({d ^ 0,i), Fig. [^a). The 
latter characteristics shows the dynamics of the contact 
line. We see that the evolution of the system consists 
of two interchanging regimes. During the time intervals 
characterized by supercritical values of the contact angle, 
\j{t)\ > 7o, the contact line keeps sliding over the sub- 
strate. This motion takes place until 7 enters the sub- 
critical domain, —70 < 7(t) < 70, where the contact line 
becomes "frozen." In Fig. ©a), the time intervals of the 
contact line being fixed are presented as the gray-filled ar- 
eas. As we clearly see, Co{t) = const here, whereas other 
characteristics are changing. The contact line remains 
fixed until the contact angle is outside the subcritical 
domain. After that, the contact line proceeds to move 
again, etc. Thus, the contact line dynamics corresponds 
to the periodic sliding interrupted by the intervals of be- 
ing completely frozen, or, in other words, to stick-slip 
motion. 

We now keep the value of the wetting parameter fixed, 
A = 1, and analyze how the amplitudes of Cp and Co 
change while varying the external frequency Q and the 
critical contact angle, 70. The corresponding dependen- 
cies are depicted in Fig. [31 As we see from the form of 
these response characteristics, the system demonstrates 
well pronounced resonances. Because the contact line 
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FIG. 4: The diagram of contact line motion on the plane 
7o). The solid lines are defined by the condition T{fl) — 70 
and separate the domains of oscillations with the fixed contact 
line (r < 70, in gray) and with the contact line moving in the 
stick-slip regime (F > 70). Courtesy of S. Shklyaev. 



FIG. 3: Amplitude-frequency response at A = 1 and different 
values of 70. Panel (a): Surface deviation at the pole, C^p. 
Panel (b): Surface deviation on the substrate, Co- 



motion with a finite value of the wetting parameter A is 
dissipative, the resonant amplitudes remain bounded. In 
Fig. [3l^a) we also present the nondissipative limiting case 
of A = 0. 

We next pay attention to the dependence on 70. The 
partial case of 70 = corresponds to no hysteresis so that 
at A = 1 the contact line keeps moving all the time. With 
the increase of 70, the part of period with the moving 
contact line becomes less and is gradually replaced with 
the regime with the fixed contact line. At large values 
of 70, the contact line remains fixed for the most part 
of period, which becomes equivalent to the case of small 
A. In other words, the dynamics with the fixed contact 
line dominates. At A = 0, the contact line is pinned, 
(0 = 0, and the oscillations are no longer dissipative, 
which results in the divergence of resonant amplitudes of 
Cp. We note that the case of A = 1 is characterized by the 
resonant frequencies close to Um for all 70, where cUm are 
the eigenfrequencies at A = 0, see Eqs. (fT2|) and (fT3|) . As 
a consequence, the amplitudes of oscillations at the pole, 
^p, are significantly higher than those on the contact line, 
Co- 

Let us now point out another generic feature caused 
by the contact line hysteresis. We start with the limit- 
ing case of no hysteresis, 70 = 0. We recall that in this 
case^ the contact line remains fixed, Co = 0, at certain 
values of the driving frequency, Q = ^ar, and any value 
of the wetting parameter, A. For this reason, the val- 
ues ilar are referred to as antiresonant. Such frequencies 
are well recognized in Fig. [S^b). As becomes clear from 
the figure, the contact line hysteresis, 70 7^ 0, transforms 
the discrete number of antiresonant points into antireso- 
nant bands of finite width. Fig. [3ljb). With the growth of 
7o, the islands of stick-slip dynamics become narrower, 
whereas the regions of behavior with the completely fixed 
contact line widen. 

We emphasize that the width of the antiresonant bands 



is determined solely by the value of 70 and is indepen- 
dent of A, which is explained as follows. The dynamics at 
frequencies within the antiresonant band corresponds to 
the oscillations with the fixed contact line, as if A = 0. It 
is clear from Eq. ([5]), that here we have r{fl) < 70 with 
r = max7(t), whereas outside the domain of antireso- 
nant behavior the opposite equality holds, r(ri) > 70. 
Hence, the border between the domains of stick-slip dy- 
namics and behavior with the fixed contact line is de- 
fined by the equality T(fl) = 70 or even much simpler: 
ro(il) = 7o, where Tq = r|A=o- As we see, the ques- 
tion about the width of the antiresonant bands can be 
efficiently answered within the nonhysteretic model.^ A 
corresponding diagram is shown in Fig. 3) 

The next question to answer is if one might expect 
any significant difference in the dynamics for other values 
of the wetting parameter. The case of small values of 
A is out of interest because the dynamics becomes very 
similar to the case of the pinned contact line, A — > 0. Our 
numerical tests show that this limit is practically reached 
at A = 1/2. As a result, the variation of 70 demonstrates 
almost no significant changes and hence the case of A < 1 
brings basically nothing new. 

Much more promising is the opposite situation, A 3> 
1. In the nonhysteretic (70 7^ 0) case, the contact line 
is not fixed and the sliding motion is predominant. In 
contrast to the case of A <C 1, interaction of the drop 
with the substrate is weakened. The eigenfrequencies 
become close to the eigenfrequencies of the even modes 
for a spherical drop, In the case of hysteresis, 70 7^ 0, 
the system is switched between two weakly dissipative 
kinds of oscillations. As we have seen for the case of 
A = 1, the stage of evolution with the fixed contact line 
is characterized by the resonant frequencies ujm- For the 
stage of sliding contact line, the resonances are found 
at Qn < Thus, a competition of the qualitatively 
different resonances is expected for the stick-slip motion. 

Our computations indicate that the described scenario 
with A ^ 1 can be observed already at A = 5. The cor- 
responding response characteristics Cp and Co are shown 
in Fig. O where the competition of pairs of neighboring 
resonances for Cp is well seen. As in the case of A = 1, the 



8 



12 
9 

a 6 

J' 

3 



4 

3 



(a) ; 




; 












A 








/ \ 








/ \ 5 
/ \^ 




// .V N 


















A = 





FIG. 6: Characteristics of the steady-state oscillations at 
A = 5, 7o = 3, and f2 = 9.8. Panel (a): Evolution of the 
contact angle, 7(t). The horizontal dashed lines are the lines 
7 = ±70. The filled areas display the time intervals of the 
subcritical contact angles, I7I < 70. Panel (b): The Fourier 
power spectrum evaluated for 7(t). 



FIG. 5: Amplitude-frequency response at A = 5 and different 
values of 70. Panel (a): Surface deviation at the pole, C^p. 
Panel (b): Surface deviation on the substrate, Co. 



growth of 7o demonstrates convergence to the resonant 
curve corresponding to the fixed contact line, A = 0. The 
transition is however nontrivial. In contrast to Fig. [3ja), 
the curve associated with 70 = has resonant peaks at 
il = ri„. As the parameter 70 is increased, the peaks do 
not simply shift from 17 = f2„ toward the values f2 = w„. 
This transition is accompanied by the emergence of in- 
termediate local maxima. Another distinction is that 
the amplitudes of resonant peaks change now nonmono- 
tonically with the increase of 70. For instance, consider 
the amplitude of in a vicinity of the first resonance, 
f2 G (2,6). At 70 = 0, we have the maximum value 

~ 7.30. With the growth of 70, the resonant ampli- 
tude starts to decrease and approaches its minimal value 

~ 3.09 at 7o = 1.45. The further increase of 70 leads 
to the growth and then divergence of the resonant am- 
plitude, as in the case of A = 0. 

The dependence of the amplitude Co on 70 and is 
qualitatively the same as described for A = 1 and is 
in agreement with the diagram of contact line motion. 
Fig. m The amplitudes of Co have, however, higher 
values because of weaker dissipation than those at A = 1 . 

We next examine the evolution of 7 and its Fourier 
spectrum evaluated at A = 5 and 70 = 3, see Fig. [HI As 
we see, the dependence 7(t) looks not only more com- 
plicated, but qualitatively different in comparison with 
that given in Fig. [21 To avoid any confusion, we here- 
after stick to the following convention. We focus on half 
the period of the signal ^(t) such that 7 > 0. Note that 
in the case discussed in Fig. [21 a) we see a single maxi- 
mum. If we now go back to Fig. [H^a) we detect the birth 
of the second local maximum, the origin of which is dis- 
cussed in a few lines. As a result, the power spectrum 
becomes wider than in Fig. [2lb) and the contribution of 
higher harmonics is stronger. It is important to indicate 
that a very similar feature has been recently observed 



experimentally, see Fig. 11, Ref. 0. Along with the ex- 
perimental study, the authors have suggested a simple 
theoretical model. Although their model is able to qual- 
itatively explain the existence of the stick-slip motion, it 
fails to reproduce the two-maxima feature in the evolu- 
tion of the contact angle, 7(t). Although the considered 
problems are not exactly the same, the advantage of our 
approach comes into play. Our model allows us not only 
to describe the stick-slip motion itself but also to capture 
the subtle feature of non-single maximum in the evolu- 
tion of contact angle. 

To get a deeper insight into the two maxima phe- 
nomenon, we provide Fig. [7| evaluated at A = 5, 70 = 3, 
and il — 11.4. We now consider half the period of Cp(i) 
with Cp > 0. Figure [TJ^b) additionally presents the pro- 
files of the drop at different times as indicated with circles 
in Fig. [Ha). 

As we see in Fig. [71[a), the dependence C,p(t) possesses 
two local maxima. One maximum, which is similarly 
present for Q-pit) in Fig. [21 concerns the stage of oscilla- 
tions with the fixed contact line. The second maximum 
is new, it relates to the stage of sliding contact line. Note 
that each stage of motion is characterized by its own res- 
onance, which are competing as we discussed for the case 
of A = 5, Fig. [5l At 7o = 0, the motion with the fixed 
contact line is characterized by the resonances at the fre- 
quencies LOrai whereas for the slip motion the resonances 
at il„ become important. At = 11.4 and 70 = 0, the 
closest eigenfrequencies are ^2 — 8.49 (fixed contact an- 
gle) and UJ2 = 10.6 (fixed contact line). If we look how 
these resonant values change as 70 is increased, we find 
that at 7o = 3 the resonances take place at ^2 ~ 9.10 and 
0)2 ~ 12.0 and the value SI — 11.4 is well in between and 
close to both of them. This explanation may also reveal 
the reason behind the two maxima as in Fig. [6l Despite 
both those maxima are found during the slip motion, 
they are caused by the competing resonances of different 
nature, as we described. 
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FIG. 7: Characteristics of the steady-state oscillations at A = 
5, 70 = 3, and Q. = 11.4. Panel (a): Evolutions of the free 
surface deviation at the pole, C,p{t) and the contact angle, 
7(t). The horizontal dashed lines are the lines 7 = ±70. The 
filled areas display the time intervals of the subcritical contact 
angles, I7I < 70. Panel (b): Profiles of the free surface shown 
at the consecutive moments of time as indicated by points 1, 
2, 3, and 4 in panel (a). 

V. CONCLUSIONS 

We have considered the dynamics of an oscillated ses- 
sile drop of incompressible liquid and focused on the con- 
tact line hysteresis. The solid substrate is subject to 
transverse oscillations, which are assumed small ampli- 
tude and high frequency. We admit that the drop is 
so small that its shape is not distorted by gravity and 
hence the equilibrium drop surface is hemispherical and 
the equilibrium contact angle equals 7r/2. To take into 
consideration the contact line hysteresis, the boundary 
condition suggested by L. M. Hocking is applied, see 
Eq. ([1]). This boundary condition involves an ambigu- 
ous dependence of the contact angle on the contact line 
velocity. More precisely, the contact line starts to move 
only when the deviation of the contact angle exceeds a 
certain critical value. As a result, the stick-slip dynam- 
ics can be observed: the system is periodically switched 
between the states with the sliding and the completely 
fixed contact line. 

The solution of the boundary value problem is pre- 
sented as series in the Legendre polynomials. Techni- 
cally, the problem is treated by building two separate 
solutions valid at subcritical and supercritical values of 
the contact angle. These solutions are different and cor- 
respond to oscillations with the completely fixed and the 
moving contact line, respectively. For the fixed contact 
line, the problem admits an analytical solution obtained 
earlier)^ In the situation with the moving contact line, a 
set of ordinary differential equations is obtained for ex- 
pansion coefficients, which are integrated numerically. At 
the critical values of the contact angle, the matching of 
the different solutions is performed. This procedure al- 
lows one to obtain the solution of the formulated problem 
at any moment of time. 

Because of dissipative nature of the Hocking condition, 
the regime with steady nonlinear oscillations is reached. 
We have measured the deviations of the free surface on 



the substrate and analyzed the frequency response at dif- 
ferent values of the wetting parameter. A, and the critical 
contact angle, 70. It is known that in the nonhysteretic 
limit, 7o = 0, no contact line motion exists at certain fre- 
quencies VL — VLar^ which are independent of the wetting 
parameter. For this reason, the values fiar are referred to 
as antiresonant frequencies. We have shown that the con- 
tact line hysteresis, when 70 7^ 0, transforms this discrete 
number of into antiresonant frequency bands of finite 
width. With the growth of 70, the parameter domains 
of the stick-slip dynamics become narrower, whereas the 
one with the completely fixed contact line grows. 

We have analyzed similar frequency response for the 
deviation of the free surface at the pole of the drop. Here, 
at relatively small values of the wetting parameter. A, res- 
onant amplification of oscillations is found at frequencies 
Q. ~ LOn for all 70, where w„ are the eigenfrequencies of 
the problem with the pinned contact line, A = 0. 

At higher values of A, the interaction with the sub- 
strate is weakened. In the case of no hysteresis, 70 = 0, 
the eigenfrequencies are close to the eigenfrequencies 
of the even modes for a spherical drop, r2„. We have 
demonstrated that the contact line hysteresis leads to 
a nontrivial shift of resonant frequencies from ~ r2„ 
to ujn as 7o — > 00. For moderate values 70 — 0(1), 
the switching between two weakly dissipative kinds 
of oscillations: with the sliding and the completely 
fixed contact line. These stages of stick-slip motion are 
characterized by the resonant frequencies r2„ and 
respectively. As a result, in the interval of frequencies 
f2 S a competition of the two resonances occurs 

and nontrivial effects can be found. Particularly, the 
evolution of contact angle has displayed the emergence 
of an additional local maximum at half a period, which 
is reminiscent of recent experimental observations^^ see 
Fig. 6(a). 
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